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SUMMARY 


The problem of chordwise, or camber, divergence at transonic and 
supersonic speeds is treated with primary emphasis on slender delta 
wings having a cantilever support at the trailing edge. Experimental 
and analytical results are presented for four wing models having apex 
half-angles of 5°, 10°, 15°, and 20°. A Mach number range from 0.8 
to 7-5 is covered. 

The analytical results include calculations based on small-aspect- 
ratio theory, lifting-surface theory, and strip theory. A closed-form 
solution of the equilibrium equation is given, which is based on low- 
aspect-ratio theory but which applies only to certain stiffness dis- 
tributions . Also presented is an iterative procedure for use with 
other aerodynamic theories and with arbitrary stiffness distribution. 


INTRODUCTION 


The current trend toward the use of thin, low— aspect— ratio, all- 
movable control surfaces on aircraft and missiles has introduced the 
possibility that divergence rather than flutter may be the primary 
aeroelastic problem for such surfaces. This possibility results from 
the fact that on all-movable surfaces the forward portion of the surface 
may be supported from the rear rather than from the side. The resulting 
chordwise divergence is characterized by a camber type of deformation 
rather than a twisting or torsion of the wing span. The present paper 
considers the problem of chordwise divergence at transonic and super- 
s °nic speeds and is principally concerned with the divergence of slender 
delta wings having a cantilever support at the trailing edge. 

The problem of chordwise bending in two-dimensional supersonic 
flow has been treated by Biot (refs. 1 and 2) and by Miles, according 
to reference 3- The three-dimensional case has also been examined by 
Miles (ref. 3) and in the more recent investigations of references 4 
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and 5- The chordwise divergence of an all-movable control at transonic 
and supersonic speeds has been treated in reference 6. 

The purpose of the present paper is to give the results of a 
reexamination and extension of the material presented in reference 5 
which had only limited distribution. The authors of this paper, who 
collaborated in preparing the results of reference have obtained 
additional results on the basis of lifting-surface theory for super- 
sonic flow. The investigation includes the development of an iterative 
solution to the equilibrium equation applicable to divergence studies 
of configurations with arbitrary stiffness distribution and, also, a 
closed-form solution applicable to certain special distributions of 
stiffness. Experimental results are presented for a series of canti- 
lever delta wings having apex half-angles of 5°> 10°, 15 °, and 20° over 
a Mach number range from 0.8 to 7-3* Comparative divergence calcula- 
tions for these cases based on small-aspect-ratio theory, lifting- 
surface theory, and strip theory are made. 


A ij 
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SYMBOLS 

curvature influence coefficient, (ft-lb)”^ (eqs. (30)) 

weighting factor in series form of Ap(x,y) (eq. (3)) 

slope influence coefficient, lb“^, (eqs. (37)) 

matrix defined by equation (42), ft^ 

matrix defined by equation (46), ft 

torsional spring constant, ft-lb /radian 

chord at wing plane of symmetry, ft 

differentiating matrix, ft“^ (see eq. (44)) 
modulus of elasticity, lb/sq ft 
local aerodynamic force per unit span, lb /ft 
concentrated aerodynamic load at 

local vertical displacement, ft 


*A> lb 



Lf\ 00 0J 
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I(x) section moment of inertia, ft^ 

integrating matrix, ft (see eq. (32)) 

J v , Jji Bessel function of the first kind of order v and p, 
respectively 

K, K(M,k,x-x* ,y-y ! ) kernel function of integral equation, ft~^ 

K q section moment of inertia constant (eq. ( 16 ) ) 

K]_, K 2 divergence constants, defined by equations (25) and (48), 
respectively 

k reduced frequency 

kj_ weight constant (see eq. ( 52 )) 

Lnm* L nm( x >y) lift function in series form of Ap(x,y) (eq. ( 5 )) 

M Mach number 

Mg(x) local bending moment, ft-lb 

q dynamic pressure, pV^/ 2 , lb/sq ft 

dynamic pressure at divergence, lb/sq ft 

Sj area of jth wing segment, sq ft 

s local wing half- span, x tan e, ft 

t local thickness, ft 

tp thickness at midspan of trailing edge, ft 

V velocity, fps 

W wing weight, lb 

x,y,z,x f ,y T Cartesian coordinates, ft (see fig. 1) 

x^,Xj chordwise coordinates of centers of ith and jth segments, ft 
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*P = c 0*i 


p = yM^ - 1 

Ap(x,y), Ap(x',y') local pressure difference, lb/sq ft 

Ax chord of wing segment, ft 

e apex half-angle (see fig. 1) 

0 angular spanwise variable defined by equation (l4) 

A characteristic parameter (see eq. (20)) 

U = 3/(3-n) 

v = n/(3-n) 

| dimensionless chordwise variable, x/cq 

dimensionless chordwise coordinate of midpoint of jth 
segment, xj/Cq 

f dimensionless chordwise coordinate of wing support station 

S P 

II T „ definite integrals in the lifting-surface-theory problem 

(see eq. (4)) 

f i ) o 

Ilnui integrals defined by equation (8), ft 

p fluid density, slugs /cu ft 

p material density, lb/cu ft 

Superscripts : 

n, m denote powers of the chordwise aid spanwise variables, 

respectively, in the expression for the thickness distribu- 
tion (see eq. (12)) 

Matrix notation: 


square matrix 


ro cq\j. M 



ro vo od \r* 
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L J 
{ } 

[ ] 


row matrix 
column matrix 


diagonal matrix 


FORMULATION OF TEE PROBLEM 


The Equilibrium Equation 

The problem at hand is that of the divergence of a delta-wing plan- 
form, shown In figure 1, which is restrained along a section normal to 
its root chord at some distance Xp from its apex. The wing is con- 
sidered to be capable only of chordwise bending - that is, bending of 
the mean camber line. The effect of spanwise variation of bending would 
be expected to be small for the narrow models with the thickness distri- 
butions studied herein. The wing in its neutral position is assumed to 
lie nearly in the xy-plane of an x,y,z coordinate system which moves 
with the planform in the negative x-direction at uniform velocity V. 

It is assumed that the wing, so restrained, will obey the mechanics of 
a simple beam and that the only external forces acting upon it are aero- 
dynamic forces which arise when it is perturbed from its neutral posi- 
tion. Under these conditions the state of equilibrium may be expressed 
by the following differential equation: 


dx 


E(x)l(x) 


dx^J 


= F 


( 1 ) 


where E(x) is the modulus of elasticity, l(x) is the section moment 
of inertia, h is the vertical displacement of any section from its 
neutral position, and F(x) is the aerodynamic force acting at 
station x. 

The present study includes several stiffness or thickness distribu- 
tions, leading to different forms of l(x), and considers various linear 
aerodynamic theories, from which are obtained different forms of F(x). 
Some examples are treated in which the thickness distribution is 
described by simple analytic functions of the chordwise coordinate, but 
consideration is also given to the possibility of an arbitrary thickness 
distribution. 
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When certain analytical expressions are chosen to represent the 
thickness distribution and are used in combination with the aerodynamic 
forces given by small-aspect-ratio theory or strip theory, closed-form 
solutions of the equilibrium equation can be obtained. For other combi- 
nations of thickness distribution and aerodynamic forces an iterative 
procedure is developed. Before discussing the methods of solution, the 
forms of F(x) and l(x) to be employed are considered. 


The Aerodynamic Force F(x) 

The forms of F(x) to be employed depend on the Mach number range 
to be dealt with and on the apex angle of the structure. If the wing is 
very slender and if the Mach lines lie well ahead of the leading edge 
(3 tan e « 1), small-aspect-ratio theory may be used. For Mach numbers 
up to that for which the leading edge becomes supersonic, lifting-surface 
theory based on kernel- function procedures is available. If the leading 
edge is supersonic and if spanwise variations of deflections are neg- 
lected, the loading based on linear theory is given exactly by strip 
theory, which for sufficiently high Mach numbers corresponds to first- 
order piston theory. 

In the present calculations, small-aspect-ratio theory and lifting- 
surface theory are used at subsonic speeds, and small-aspect-ratio 
theory, lifting- surface theory, and strip theory are used at supersonic 
speeds . Both small-aspect-ratio theory and strip theory are used for 
comparison throughout the range of supersonic Mach numbers considered 
even though the Mach number range in which the theories are valid is 
exceeded. 

Small-aspect-ratio theory . - An expression for the aerodynamic force 
given by small-aspect-ratio theory for the case of steady flow can be 
obtained directly from the work of Jones (ref. 7) or from Garrick's 
results for unsteady flow (ref. 8). The force can be expressed in the 
following form: 



( 2 ) 


CO C\ C\J 
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where 


q = |pv 2 

Lifting-surface theory .- In order to obtain an expression for the 
force F(x) based on lifting-surface theory, it is necessary to deter- 
mine first the chordwise and spanwise distribution of pressure acting 
on the surface. For this purpose use is made of an integral equation 
which, in steady flow, relates the slope of the wing surface to the 
pressure distribution. 

The integral equation can be written as 


dh 

dx 


1 

4jrpV 2 



Ap(x* ,y ! )K(M,k,x-x* ,y-y r )dx* 


dy ! 


where K(M,k,x-x f ,y-y j ) is the kernel function. The quantity K/4jtpV^ 
is the mathematical expression for the downwash induced at any point 
x,y by a unit force acting at any other point x*,y T . The area A over 
which the integration extends is the portion of the wing in which a 
pressure pulse must occur in order to induce vertical velocity at the 
specified point x,y. In subsonic flow A corresponds to the entire wing 
surface; in supersonic flow A is that portion of the surface bounded 
by the planform edges and the forward Mach cone from the point x,y. A 
numerical method of solving the integral equation for the subsonic case 
is described in reference 9j the procedure for the supersonic case is 
based on similar techniques but involves differences in the form of the 
pressure distribution and in the kernel of the integral equation. 

An approximate solution to the integral equation involves expressing 
the unknown pressure distribution as a sum of chosen modes of lift func- 
tions (of forms appropriate to the planform and Mach number range 

under consideration), each weighted by a constant coefficient to be 
determined. The following expression is employed: 


Ap(x,y) = 


(3) 


Through the use of this expression the integral equation can be repre- 
sented as a summation of definite integrals and can be solved by numeri- 
cal methods. The following matrix form of the integral equation can 
then be written: 


{e} - M <>} 


( 4 ) 
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where the elements of 11-^ represent the definite surface integrals 
of the products of the functions and the kernel function K and 

dir 

the elements — represent the wing slope at a number of selected 
dx 

points on the wing surface. Equation (4) may therefore be regarded as 
a set of simultaneous equations from which the values of the weighting 
factors a^ can be obtained once the definite integrals have been 
evaluated. (The evaluation of these integrals constitutes the major 
task in the kernel-function procedure and is accomplished by use of the 
methods and the computing program described in reference 9; it is not 
pertinent to the present discussion, however, and need not be considered 
herein. ) 


A premultiplication of each side of equation (4) by the 


matrix 



gives the expression 


inverse 



(5) 


which can be used to obtain values of the weighting factors for 

any prescribed slope distribution. Once these weighting factors a M 
have been determined, the pressure Ap(x,y) associated with the pre- 
scribed slope distribution is defined by equation (3) . 

For use in the influence-coefficient procedures to be described 
subsequently, it will be necessary to obtain the forces acting on each 
of the several segments (fig. 2) into which the wing is divided; this 
will require an integration of the pressure distribution over each 
segment. The force on the jth segment Fj may be expressed as 


(6) 

(7) 


II ™ ) = ff s L nm (x »y) dx d y 

and where Sj denotes the area of the jth segment. 




Ap(x,y)dx dy 




or, applying equation (3), 


Fj = 8rrq 


II 


0) 

'run 


{v) 


where 


(8) 


^ mao oj 



Str ip theory,- If the leading edge is supersonic and if there is 
no spanwi.se variation of wing deflection h, the aerodynamic loading 
is given exactly (within the limitations of linear theory) by strip 
theory, (See, for example, refs. 10 and 11.) The force F(x) for 
use in equation 0 ) is expressed, as 


or 


F(x) 


2 P v 2 r x tan e dh dy 

P J -x tan e dx Y 


F(x) = - §3= x tan e — (9) 

' P dx 

If M » 1, F(x) can be expressed approximately as 

F(x) = - x tan e — (10) 

' M dx 

which is the result given by piston theory (ref. 12) for the case of 
zero th i ckne s s . 


Area Moment of Inertia l(x) 

The term l(x) in equation (l) is the section area moment of 
inertia and is therefore directly related to the thickness by the 
expression 

p x tan € nt/2 p x tan e 

l(x) = / dy / z 2 dz = — / t^dy (ll) 

' J -x tan e J -t/2 12 J -x tan e 


where £ is the apex half -angle. 

As noted previously, it has been found that, when certain forms 
of the thickness distribution are used in combination with aerodynamic 
forces based on small-aspect-ratio theory or strip theory, closed-form 
solutions of the equilibrium equation are possible. A fairly general 
form of the thickness distribution which leads to an exact solution is 
the following expression: 
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t 



( 12 ) 


where t Q is the thickness at the point x=cg,y=0 and where m and 

n are positive integers. This expression for t contains a chord- 

X. 

wise variation of thickness given by the term in — and a spanwise 

c 0 

variation in thickness given by the term involving ^ . It should be 

O 

remarked that for n > 3 the thickness distribution is characterized 
by a cusp at the apex so that at the apex the leading edge is infinitely 
sharp. Such thickness distributions are not considered herein since 
linear beam theory would not be expected to apply. 


In order to proceed toward an expression for the moment of inertia 
l(x), equation (12) may be substituted into equation (ll) to obtain 


I(x) 



(13) 


where s = x tan e is the local wing half-span and where symmetry 
about the center line is assumed. The integral in equation (13) can 
be expressed conveniently in terms of the gamma function and for this 
purpose' there is introduced the angular variable 6 defined as 

0 = sin" 1 ^ (14) 

s 

in terms of which equation ( 13 ) becomes 


l(x) 



/" 

Jo 


*/2 (2+m) /2 

(cor. 0) d0 


(15) 


If the integral in equation (15) is denoted by Kg, it may be shown 
that 


Kq - 


>*/2 ( 2 + m )/2 

(cos 9) d0 = 


0 



(1 6 ) 


£ 


ro cx> vji f 



hi Lf\ CO (\J 
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where T denotes the gamma function. The inertia term l(x) may 
then be expressed as 



( 17 ) 


Numerical values of the section moment-of-inertia constant Kq for 

several values of m are given in table I. Included with table I are 
sketches of the sections to which the various values of m pertain. 


CLOSED-FORM SOLUTIONS OF THE EQUILIBRIUM EQUATION 


Exact solutions of the equilibrium equation can be obtained for 
a few special choices of the section area moment of inertia l(x) 
given by equation ( 17 ) and for the aerodynamic -force distributions 
F(x) given by small-aspect-ratio theory or by strip theory. These 
solutions can be used to assess the accuracy of the more generally 
applicable iterative procedures to be developed. Before considering 
these solutions it is convenient to arrive at a nondimensional form of 
the equilibrium equation by introducing the change of variable 

x = c Q t 

so that equation ( 17 ) becomes 



With this expression for I^cq£), equation (l) may be written in the 
desired nondimensional form as 


VM 

"6\ c o/ 


3 


tan e 






( 18 ) 


If the expression for F^c^l) given by small-aspect-ratio theory 
(eq. (2)) is employed and if E(c 0 i) is taken to be constant, equa- 
tion (l8) can be written as 



12 


o 

d.- 




+ 


A 2 




= 0 


where 


(19) 


^2 _ 12-.it q tan £ 



( 20 ) 


Solutions to equation (19) will be obtained fcr values of n of 0, 1, 
2, and 3 for two types of wing mounting. 


First, however, it is noted that a solution of equation (l8) can 
also be obtained, provided n = when the aerodynamic -force distribu 
tion is given by strip theory (eq. (9)). In this event equation (l8) 
becomes 


L 


o 


K of t o 
6 V c o 



e 


dr 


■ it 


(Mb 


4 d^h | 
d| 2 j 


8q dh 

+ | tan e — 

0 d£ 


or taking E(cq§) as a constant and performing the indicated differen- 
tiations of the first term results in 


F + + 12| 2 + At 

4 , t 3 * e 2 d£ 


= o 


dr 


dr 


dr 


in which 



This form of the equilibrium equation is of the form of Cauchy's equi- 
dimensional linear differential equation, the solution of which is 
discussed, for example, in section 1.6 of reference 13 . 


LT\ CO 



Co '• 


L ' 1 ^ ^ 00 ^ion 1 t < ii e 1 l r a. j 1 i ng - • E dg e - C ai j. t i le ve r e d D e 1 la 

wiiig Baaed or:- Small-Aspect -Ratio Theory 

la uidci’ lu complete the formulation ol l the divergence problem as 
f. oS^d by equation (l-9)j it is necessary to impose certain condi- 
tions at the apex and at the point of support. Conditions to be imposed 
. at the a P ex irn Piy no bending moment or shear and are 


li .111 £ 


n+I d 2 h 

d^ 2 




. 2 , 


lim 


d ~h \ 


-t-J-ui c 

i ->0 \ d ^2 


\ ° A ~ L " in " tx ailing-edge-ca.nti levered. delta wings, the conditions to be 
imposed imply that there is no displacement or slope at the built in 
edge and. are 


fool 


/dh\ 

UiA 


(22) 


&=1 


L^r one integration and the imposition of the apex conditions 
eqnat.-i ,.>n ( ! 9 ) becomes 


o 

dlh 


\ 


n+ 1 sLl aV li = 


die 


d£ 


0 


t r if dh , 

’ ' d| ' 


d 4ml df \ 
d|\ dp 


+ A^pf - 0 


’Thic may be recognized as a special form of Bessel 
for examjrle, eq. (l23e), p. 167 , of ref. Id.) If 
under the ep&x conditions imposed by equation (21) 


equation. (See, 
r 5 > a so Jet-ion 
found to be 
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h = C n 


r n/2 J v 


2> ,(3-n)/2| 


n 


t 


d£ 4- C 2 


( 25 ) 


where J v is the Bessel function of the first kind of order v 
(v = n/(3 - n)) and and C 2 are constants. (if n = 3, the 

problem is soluble as an equidimensional equation (sec. 1.6, ref. 13 ) • 
From equations (22) and (23) , the eigenvalues for the cantilevered case 
(n / 3) are the roots of 


27\ 


3 - n 


(ft) 


(3-n)/2 


= 0 


(24) 


In order to assess the accuracy of the iterative procedures to be 
discussed subsequently, solutions to the differential equations have 
been obtained for constant E and for a thickness distribution defined 
by setting m = 0 in equation (12) so that 


t 


- t o( £) 


n/3 


and for n = 0, 1, 2, and 3- For n = 0, 1, and 2 the roots of the 
Bessel functions in equation (24) and the corresponding values of ?\ 
are as follows: 


n 

n 

Equation 
in A 

Root of 
J v = 0 

A 

' 3 - n 

0 

0 

J °(f) ■ 0 

2.4048 

3.6072 

1 

1/2 

Ji/ 2 (A) = 0 

It 

rt 

2 

2 

J 2 (2A) = 0 

5.1356 

2.5678 


For n = 3 there is found a corresponding value of A of 1.499* The 
quantity ?\ 2 in equation ( 19 ) may be used to form a convenient diver- 
gence constant K-^ defined as 


Kn 


_2 q^ tan € 

A d 



Values of for the foregoing examples are listed in table II. 


(25) 
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Exact Solution for the All-Movable Control 


Based on Small-Aspect-Ratio Theory 


An alternative configuration which might be treated is the all- 
movable control attached to a torque rod at | ^ where £ is the 

chordwise coordinate of the pitch axis and is referred to Cq as unit 


length. For this case the conditions to be imposed at the section of 
support would be 


where 


h (h) - ° 



r 1 

m b( 6 p) - c 0 2 / (5 - 6p) F (t)as 

j o 


or, from equation (2), 


( 26 ) 


M^lp) = -2rtq tan 2 ec 0 2 



(27) 


is the total aerodynamic moment about the pitch axis £ = | and 
where Cp represents the torsional spring constant of the torque rod. 


In order to obtain an expression for the eigenvalues for the all- 
movable control attached to a pitch spring, the expression for h 
given by equation ( 23 ) must be substituted into equation ( 27 )} after 
the resulting expression is integrated, the expression which defines 
the eigenvalues is 




(5-n)/2 


(28) 


where p = 5/(3 - n) and, as stated previously. 


V 


n/(3 - n) . 
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IMPERATIVE SOLUTIONS OF THE EQUILIBRIUM EQUATION 


The exact solutions of the equilibrium equations which were out- 
lined in the preceding section are applicable only to calculations 
based on certain aerodynamic -force representations and thickness dis- 
tributions. For other cases it is necessary to seek an approximate 
solution; in the present investigation an iterative procedure based on 
influence coefficients is employed. 

No proof of convergence of the iterative procedure is attempted. 
In the calculations described in reference 5 for the case of an all- 
movable control with the hinge axis not at the trailing edge, some 
difficulty with convergence was experienced in a few cases for which 
the control surface was very stiff and the torque rod was very weak. 

It was found, however, that by averaging the results of two successive 
cycles of iteration and employing the average to start the next cycle 
convergence could be obtained. No such difficulty was encountered in 
the present calculations for the delta wings with a cantilever support 
at the trailing edge. 

Formulations of the equilibrium equati based both on curvature 
influence coefficients (eq. ( 31 )) and on sl^pe influence coefficients 
(eq. (34)) will be given. Either formulation may be used, the choice 
being dictated by convenience In obtaining :he necessary structural 
and aerodynamic data , 


Influence-Coefficient Forms of the Equilibrium Equation 

As is the usual practice , the wing is considered to be made up of 
a finite number of segments (fig. 2). In tie calculations for super- 
sonic flow, 10 segments were used and comparisons with exact results 
(to he discussed in a subsequent section) slowed this number of seg- 
ments to be adequate. In the calculations "or subsonic flow, first 
10 and then l6 segments were used; the additional segments were added 
near the apex between 0 and 0.3 chord. The results were not changed 
by the addition of these segments. 

The force acting on each segment is assumed to be concentrated at 
the midpoint of the segment and the slope o: ' curvature is measured 
there. The analysis proceeds from the fundamental beam relation between 
curvature and bending moment, namely, 



Mg(x) 

El(x) 


where M-g(x) is bending moment. 


( 29 ) 



H-l \r\co cm 


3X 
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Curvature equation .- For arbitrary distributions of stiffness, 
equation ( 29 ) can be used directly to formulate an iterative procedure 
based on curvature influence coefficients. If denotes the curva- 

ture produced at Xj_ by a unit load Fj at Xj, there can be obtained 
from equation ( 29 ) the relations 


dx2 L ± 

A. . = — 

1J F j 


x i - X J 

E( x i)I( x i) 


( x i > Xj) (30a) 


and 


A- 


ij 


= 0 


( x i ^ Xj) (30b) 


Equations (30) are used to evaluate the curvature influence coefficients 
A±j and also to provide the following matrix form of the equilibrium 
equation: 



The aerodynamic forces Fj acting on the various segments can be 

expressed as functions of the curvature to obtain a form of the equilib- 
rium equation which may be iterated until convergence upon a curvature 
distribution is obtained. The manner of doing this for the various 
aerodynamic theories is indicated in the next section. 

Slope equation .- A matrix form of the equilibrium equation based 
on slope influence coefficients may be obtained directly from equa- 
tion ( 31 ) by the introduction of a suitable integrating matrix pjJ 

which gives the results 
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where an element B i j denotes the slope produced at by a unit 

load at x j . 

A sample integrating matrix appropriate to the 10-segment system 
and the trailing-edge-cantilevered models treated herein is given in 
table III. Multiplying both sides of equation (31) by the matrix ^ 

leads to the following expression: 



NhlW ■ N( f j) 


( 54 ) 


This expression provides the basis for iteration on the slope distribu- 
tion for a structure with arbitrary stiffness. 

For the special cases for which the thickness distribution is 
specified by equation (12) and for which the material is homogeneous, 
equation (29) can be integrated in closed form to obtain an analytic 
expression for the slope influence coefficients j • For such cases 

l(x) is given by equation (17) in terms of the dimensionless variable 
Z = x/c 0 as 

I (x ) = K ° 6 - °- (|) n c 0 l tan e (35) 


By substituting this result into equation (3Ca.) the following expres- 
sion is obtained for the curvature at | due to an arbitrary load F* 

at where | > |j: 


A . 6f j s - h 

dx 2 KQEt 0 5 tan e (l/ n+1 


This result can be integrated to obtain the following expressions for 
B^j for n = 0, 1, 2, and 3 and for Z± ^ t ^ - 


ro oo vji t-* 



I\J UJ\J\ 
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For n = 0, 


dh\ 

dxi 


B U 


6c 


0 


' j KQEtQ^tan 


(1 - 6 ± + Sj log £ ± ) (57a) 


n = 1, 


B 


ij 


6c o /. _ , . _ h} 

K^EtQ^tan € \ ^ 1 Si/ 


(57b) 


n = 2, 


B 


6cr 


ij 


' Sj i fij 

-l + + -±- - J 


K^EtQ^tan € \ 2 ^ 


(57c) 


and n = 3, 


6c, 


B ij - 


Ll 


KQE^tan € ^ 2 5 2l ± 2 3i ± 3 


(57d) 


For | . < i . , 


B ij " B jj 


(58) 


Solutions Based on Small-Aspect-Ratio Theory 

Equation (2) can he used to express the aerodynamic -force distribu- 
tion given by small-aspect-ratio theory in terms of slope and curvature 
as 


F(x) 


-2nq tan e 


2x^ + 

V dx 



(59) 
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If the force is evaluated at Xj, the center of the jth segment, and is 
considered to be constant over the segment, the matrix of aerodynamic 
forces <^F for use in the equilibrium eqvation may be then written 
as 



where Ax is the chord of the (equal) segments. 

Curvature equation .- For use in equation (31) it is desired that 
Fj be expressed as a function of curvature only. For this purpose 

equations ( 32 ) and (40) can be used to write 


where 




Substitution of equation (4l) into equation (31) yields 


(41) 


(42) 



which may be solved by iteration upon the curvature to obtain the value 
of dynamic pressure at divergence . 

Slope equation . - It is possible to obtain the matrix of 

equation (kO) in a form which involves only the slope and which can 
therefore be used with slope influence coefficients in equation (^4) . 

For this purpose a differentiating matrix [d] is needed to yield 
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A suitable differentiating matrix for use with the 10-point system of 
the present study was obtained in reference 6 and is given in table IV. 
The use of this differentiating matrix in equation (40) gives 



( 45 ) 


where 


B2I = tan^e 


MW + 2 N] 


(46) 


Equation (45) may be substituted into equation (34) to obtain the fol- 
lowing expression: 



(47) 


which may be solved by iteration upon the slope to obtain the value of 
dynamic pressure at divergence. 

Numerical results .- In order to help determine the accuracy of the 
iterative procedure, solutions to equations (43) and (47) have been 
obtained for the thickness distributions for which closed-form solutions 
were obtained. 

Values of the divergence parameter K-^, equation (25 ) , obtained 

from a 10-station representation of the wing are compared in table II 
with values from the closed-form solution and are seen to be in very 
good agreement, which indicates that 10 stations are adequate. It is 
noted that the force distribution, and consequently the divergence 
parameter, determined by low-aspect-ratio theory is independent of Mach 
number . 

One further check on the accuracy of the iterative solutions based 
on partitioning the wing into discrete segments was carried out by 
assuming a polynomial expression for the displacement that satisfied 
the end conditions and by applying the method of Stodola and Vianello 
(sec. 5 . 5 , ref. 13) to the loading equation (l 8 ) . The strip-theory 



22 


representation of the aerodynamics was used. With this representation 
of the loading, another convenient divergence constant can be 

defined as 


K 2 



m 


Iteration on the loading equation with the polynomial expression for 
the displacement led to a value of of 60.59^- Iteration on the 

slope equation for the 10-segment, partitioned wing gave a value of 
K 2 of 60.241. 


Solutions Based on Strip Theory 

The forces on the various wing segments given by strip theory, 
equation ( 9 ), may be written for use in the slope-influence-coefficient 
equation (3^) as 


{Fj} = - f Axq tan 



the substitution of which in equation (3*0 yields 


n 


i&\ 


, \<lx/ 




- — Axq tan e 

P 


NNPJ 


(*+9) 


Alternatively, equation (32) may be used to express the forces in 
terms of curvature as 



o 

— Axq tan 

p 


so that equation (31) becomes 




( 50 ) 
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Either equation (49) or equation (50) may be iterated to convergence 
to obtain the critical value of q. Results obtained for the flat-plate 
delta wings of the present study are given in a subsequent section. 


Solutions Based on Lifting-Surface Theory 

The aerodynamic forces given by lifting-surface theory have been 
used to develop an iterative solution of the equilibrium equation based 
on slope influence coefficients. It could be readily modified for use 
with curvature influence coefficients. 


It is recalled that the pressure distribution over the wing is 
given by a series expression with arbitrary weighting factors a^ 
which are defined by equation (5) for any specified slope distribution. 
In concept, one uses equations (5) and ( 7 ) to obtain the forces Fj as 


or 




( 51 a) 


( 51 b) 


and substitutes this result into equation (3*0* This provides a form 
of the equilibrium equation upon which iteration to convergence can be 
performed . 


In practice, it is usually necessary to employ several steps in 
each cycle of iteration. This is brought about by the fact that the 
inverse matrix in equation (5) is developed for a particular array of 
control points which do not, in general, coincide with the midpoints 
of the segments used in the iterative procedure. To perform one cycle 
of iteration the following procedure is used: first, assume a set of 

slopes at the control points and use these slopes in equation (5) to 
obtain a set of weighting factors next, use these weighting fac- 

tors in equation ( 51 a) to obtain the forces F y f and, finally, substi- 
tute the values of the forces into equation (34) to obtain the slopes 
at the midpoints of the segments . To begin the next cycle of itera- 
tion, slopes appropriate to the control-point locations are sorted out 
of the results of equation (34) and are used in equation (5)> and so 
forth. Usually, convergence is obtained after about four or five cycles 
of iteration. 
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DIVERGENCE EFFICIENCY OF A CANTILEVER DELTA WING 


As a matter of interest it is noted that one can determine the 
thickness distribution of the family of distributions expressed by 
equation (12) for a solid cantilever delta wing that will produce the 
highest dynamic pressure at divergence for a given total weight. 

The weight W of a solid delta wing is found by integrating the 
thickness distribution and combining the result thus found for the 
volume with the density p of the material to give an expression of 
the form 


p 

W = kjptQCQ^ tan e (52) 

where the form of the constant k^ depends on the form of the thick- 
ness distribution. This result may be used together with the closed- 
form result for the divergence parameter given by equation (25) to 
obtain the expression 


q a _ K i EK o 

p^cQ^tan^e 


Expressions of the same general form are obtained from the iterative 
solutions and from the various aerodynamic theories. 


Values of 



obtained with low-aspect-ratio theory and with 


piston theory are shown in figure 3 for the thickness distribution 
given by equation ( 12) for m = 0. The values of q^y/w^ have been 

normalized to the value for the constant-thickness wing (n = 0). The 
thickness distribution for maximum mass efficiency could be determined 
by variational procedures; however, the maximum value of q^W^ 

the associated thickness distribution would rot differ appreciably from 
the results for n = 1. 


APPARATUS AND TESTE 1 


Models 

A series of delta-wing models of constant thickness and of apex 
half-angles of 5°> 10°, 15°, and 20° were constructed. A sketch of the 
various model configurations is presented in figure 4. The thickness 
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at the base of the delta wing before attachment to the wind-tunnel sting 
was increased to simulate cantilever base boundary conditions with a 
minimum of aerodynamic interference . For the transonic tests , two 
model configurations were used. Both configurations had a 10-inch chord 
with one model constructed of 0 .048-inch-thick aluminum alloy and the 
other of 0.0385-inch-thick aluminum alloy. For the supersonic and 
hypersonic tests , wings of 0.048-inch-thick aluminum alloy with a 
6-inch chord were used. 


Wind-Tunnel Tests 

Transonic tests .- Divergence data were obtained at transonic speeds 
with the models mounted on a sting in the Langley 2-foot transonic 
aeroelasticity tunnel with Freon-12 as a test medium. The tunnel Mach 
number was held constant and the test-section density, and hence dynamic 
pressure, was slowly increased until divergence occurred. Occasional 
adjustments in the sting angle of attack were required to correct for 
changes in tunnel-flow angularity so that the model would remain at zero 
lift until the divergence dynamic pressure was reached. Small, rela- 
tively high-frequency oscillations of predominantly the apex region of 
the delta wing usually occurred at low dynamic pressure and continued 
intermittently as the dynamic pressure was increased up to the diver- 
gence condition. These oscillations were believed to be associated with 
flow separation or with flow irregularities in the tunnel stream and 
with low structural damping of the models. The divergence dynamic pres- 
sures were very sharply defined and were marked by one or two large 
excursions of the tip of the model just prior to divergence. The high- 
frequency oscillations of the tip region usually continued intermit- 
tently during the preliminary excursions and during the divergence. It 
was not determined whether these oscillations had any effect on the 
divergence characteristics of the models; however, the data faired 
fairly well with the supersonic data where the oscillations were not 
apparent. The model motion, when divergence was reached, was very 
rapid, with the deflection increasing until the model was bent past 90 ° 
to the airflow. Some of the models after testing are shown in figure 5- 

Supersonic and hypersonic tests .- The supersonic tests at M = 2.0 
and M = 3*0 were conducted in the Langley 9“ by 18-inch supersonic 
aeroelasticity tunnel with air used as a test medium. The hypersonic 
tests at M = 7*3 were made in the Langley hypersonic aeroelasticity 
tunnel which uses helium as a test medium. The test procedures at 
supersonic and hypersonic speeds were the same. The models were mounted 
on a sting and the tests were made at a fixed Mach number. The stagna- 
tion pressure was increased until the model diverged. A strain gage 
on each model was used to correlate the time of divergence with the 
recorded tunnel dynamic pressure. Each complete test lasted from 4 to 
10 seconds. There was insufficient time to adjust the model angle of 
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attack during a test and several models that were not alined properly 
slowly loaded up and failed. They were discarded and the tests were 
remade with new models . 


RESULTS AND DISCUSSION 


The transonic and hypersonic divergence data of the present inves- 
tigation are presented in figures 6 and 7- Experimental results are 
compared with results of calculations based on small-aspect-ratio theory, 
lifting-surface theory, and strip theory. A nondimens ional divergence 


parameter 



in which 




denotes the dynamic pressure at diver- 


gence, is employed in presenting the results. In figures 6(a) to 6(e), 
the divergence parameter is plotted as a function of the apex half- 
angle for the different Mach numbers. In figure 7> the divergence 
parameter is plotted as a function of Mach number for the four different 
apex half -angles used. In figures 6 and 7 the solid portions of the 
curves indicate the range of Mach numbers or apex half-angles for which 
the theories might be expected to be valid. Thus, for small-aspect- 
ratio theory the solid curves denote conditions where the Mach angle is 
at least twice the apex half-angle, and for strip theory the solid lines 
denote conditions where the component of stream velocity normal to the 
leading edge is supersonic . 


Examination of the experimentally determined values of the diver- 
gence parameter shows that for a given Mach number (figs. 6(a) to 6(e)) 
increases with decreasing apex angle and for a given apex angle 

(fig. 7) increases with increasing Mach number. With regard to 

the various analytical approaches to be considered, the figures illus- 
trate the following facts: Small-aspect-ratio theory is independent 

of Mach number and inversely proportional to tan e; strip theory is 
independent of apex angle and directly proportional to 0; and lifting- 
surface theory is a function of both Mach number and apex angle. 

At M = 0.8 (fig. 6(a)) lifting-surface theory and small-aspect- 
ratio theory give very similar results. For e = 20°, both indicate 
values of well above that found experimentally and become increas- 

ingly high relative to experiment as e decreases. Several approaches 
were taken in trying to improve the results given by subsonic lifting- 
surface theory. In obtaining the aerodynamic forces three different 
arrays of control points were used, including one distribution based 
on Gaussian techniques. In performing the iterations on the equilibrium 
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equation the l6-segment system shown in figure 2 , with close spacing 
near the leading edge, was used to examine the possibility that the 
chordwise centers of pressure might be too far removed from the centers 
of the forward segments in the 10-segment system. The results were 
essentially insensitive to these changes. 

For the ranges of Mach number and apex angle for which small- 
aspect-ratio theory might be expected to apply in supersonic flow, 
values of dynamic pressure at divergence predicted by this theory are 
generally above the experimentally determined values. When applied 
beyond the range of expected validity, small-aspect-ratio theory pre- 
dicts values of q^ well below experiment. 

For Mach numbers above about 2.0 lifting-surface theory (and strip 
theory for supersonic leading edges) predicts the experimental trends 
quite well except for the wing having an apex half-angle of 5°* For 
this case none of the theoretical approaches used gives adequate agree- 
ment, possibly because viscous effects may be quite strong. When strip 
theory is applied for subsonic leading edges, the predicted values of 
divergence dynamic pressure are generally quite conservative; first- 
order piston theory would yield values higher by the factor M/p. 


CONCLUSIONS 


The results of the analysis of the static divergence of low-aspect- 
ratio triangular wings and the comparisons of these theoretical results 
with the experimental data lead to the following conclusions: 

1. A general iteration procedure is developed for computing diver- 
gence dynamic pressures for delta wings using aerodynamic forces given 
by small-aspect-ratio theory, lifting-surface theory, and strip theory. 
The procedure can be extended to arbitrary planforms and other aero- 
dynamic repre sent at ions . 

2. Certain special thickness distributions led to closed-form 
solutions for the divergence dynamic pressure when the aerodynamic 
forces were given by small-aspect-ratio theory. The closed-form solu- 
tions were useful for evaluating the numerical procedures. 

3. At a given Mach number the experimental divergence dynamic 
pressures of the cantilever delta models of uniform thickness, which 
were tested, increased with decreasing apex angle; for a given apex 
angle the experimental divergence dynamic pressures increased with 
increasing Mach number. 
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4. In subsonic flow lifting-surface theory and small-aspect-ratio 
theory give very similar results but predict dynamic pressures at diver- 
gence which are high relative to experiment. 

5- Values of dynamic pressure at divergence given by small-aspect- 
ratio theory vary from generally above experiment in the ranges of Mach 
number and apex angle for which the theory is expected to be valid to 
well below experiment for conditions outside this range. 

6. For Mach numbers above about 2.0 lifting-surface theory (and 
strip theory for supersonic leading edges) predicts the experimental 
trends quite well except for the wing with an apex half -angle of 5° for 
which viscous effects may be important. 


Langley Research Center, 

National Aeronautics and 
Langley Field, Va., 


Space Administration, 
June 9, I960. 
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TABLE I.- SECTION MOMENT-OF- INERTIA CONSTANT Kq FOR 
VARIOUS SPANWISE SECTIONS 






TABLE II.- COMPARISON OF RESULTS 


SOLUTIONS FOR VARIOUS CHORDW] 
CONSTANT SPANWISE £ 

t(x,y) = 


n 


Qd tan 

€ 

f 

" 


Exact 

solution 

Iteration on 
curvature 
equation (31) 

Iteration on 
slope 

equation (3^) 

0 

49.5 

49.2 

49.4 

1 

37-7 

1 

37.3 

37-3 

2 

25.08 

24.49 

24.48 

3 

8.58 

8.54 

8.42 


< 


r 


L_ 

0 
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0.2541 


-0.0329 


- 0.0250 


-0.0265 


-0.0259 


-0.0263 


-0.0259 


-0.0265 


-0.0250 
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TABLE III.- INTEGRATING MATRIX 



FOR EQUATION (32) 


Common factor, 



1.9126 

-1.3974 

6.0090 

-6.6564 

9.4773 

0.6105 

0.1557 

3.8041 

-4.2750 

7.6157 

0.2446 

-0.7912 

4.4101 

-4.8225 

8.0083 

0.2671 

-1.2230 

3.6389 

-4.5240 

7.8297 

0.2602 

-1.1758 

3.1415 

-5.1803 

7.9904 

0.2642 

-1.1982 

3.2299 

-5.7500 

7 .4208 

0.2604 

-1.1788 

3.1662 

-5.5893 

6.7644 

0.2662 

-1.2072 

3.2513 

-5.7678 

7.0629 

0.2510 

-1.1555 

3.0472 

-5.3753 

6.5154 

0.3314 

-1.5058 

4.0660 

-7.2368 

8.8969 


-5-7545 4.8079 


A. 7357 4.4375 


-4.9398 4.5092 


-4.8547 4.4808 


-4.9184 4.5003 


-4.8299 4.4779 


-5.3273 4.5250 


-6.0985 4.0932 


-5A926 3-1463 


-7.6975 4.6994 


-0.4640 1.3U5 


-0.3837 1.3036 


-0.3988 1.3051 


-0.3930 1.3046 


-0.3968 1.3049 


-0.3928 1.3046 


-0.3998 1.3051 


-0.3772 1.3036 


-O.7431 1.3115 


-2.0452 


-O.0329 


1.0246 



54 


TABLE IV.- DIFFERENTIATING MATRIX [ij FOR EQUATION ( 44 ) 
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Figure 6.- Concluded* 
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